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SUMMARY 

A linearized solution is obtained for the problem of plastic 
deformation of a thin plate with a circular hole in the strain- 
hardening range. This solution is based on the deformation theory 
of plasticity for finite strains. The two final simultaneous non- 
linear differential equations obtained previously are linearized 
by replacing the nonlinear terms with a simple linear term involving 
a constant, which is determined by minimizing the effect due to the 
linearization. Stress and strain distributions of the problem for 
a given material and a given maximum strain can be obtained through 
a simple multiplication by means of the tables presented herein. 
Equations for stress and strain concentration factor are also given. 

Numerical examples are calculated by the linearized method. 
Calculations are also made for ideally plastic materials and for a 
power-law approximation. An approximate method of solution for 
stress is also presented herein. The results obtained by different 
methods are compared and the following conclusions were obtained: 

1. The results obtained by the linearized method compare closely 
with those obtained previously without linearization. 

2. The variation of a parameter determined from the octahedral- 
shear stress-strain curve of the material can be used as a simple 
general criterion of applicability of deformation theory for this 
problem. 

3. The ratios of the strains along the radius to their Tnaytmnm 
values and the ratio of principal stresses are essentially independent 
of the octahedral shear stress-strain curve of the material, but the 
stress distributions are very dependent on the material. 
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4. The solution for ideally plastic material with the infinites- 
imal strain concept gives good approximate values of strains , hut 
not of stresses. 

5. Sufficiently accurate values of principal stresses can he 
obtained hy using the strains determined by. the ideally plastic 
approximation together with the actual octahedral shear stress-strain 
relation of a given material. 

6. If a simple analytical function representing the octahedral 
shear stress-strain relation is required for analysis, the power- law 
approximation can he used. 


UTOBODUCT ION 

In the design of highly stressed machine parts, knowledge of 
the stress and strain concentration due to a hole and also of the 
distributions of stresses and strains in the strain- hardening range 
is desirable. The solution for a thin plate with a circular hole 
for ideally plastic material under infinitesimal strain has been 
obtained by Nadai (reference 1, page 189). The same problem has 
been solved (reference 2) in the strain- hardening range by using 
the deformation theory of plasticity and employing the finite- 
strain concept; numerical results were obtained for two materials 
in reference 2. Because of the nonlinearity of the equations, 
however, the results were more or less limited to the two particular 
materials used in the calculation. In order to obtain general 
information applicable to most materials, it is desirable 'to inves- . 
tigate the possibility of linearizing the equations. A linearized 
solution will also greatly reduce the amount of calculation required 
for a given case. 

Generalization of these results also has theoretical signif- 
icance in the theory of plasticity in that a general criterion of 
the applicability of the deformation theory may be obtained. This 
theory, formulated by Hencky and Nadai (reference 3), states that 
the state of stress is a unique function of the state of strain when 
the directions and the ratios of the principal stresses remain constant 
during loading in the absence of time and temperature effects and of 
unloading. It is true for plastic deformations, which are irre- 
versible processes, that the plastic strains at a certain state depend 
on the path by which that state is reached. For a particular path, 
the relations between the state of strain and stress are unique if 
the material is isotropic before loading, as is generally assumed. 
Experiments on thin tubes (references 4 to 10) show that the 
relations between the strain and the stress for the different 
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paths, each under a constant ratio of principal stresses (also con- 
stant directions), are close. The results obtained by the deformation 
theory therefore have physical significance if the previously mentioned 
conditions (the directions and the ratios of principal stresses remain 
constant during loading) are satisfied. Ho inconsistency exists in 
the deformation theory if it is used for the cases where these con- 
ditions are fulfilled. Even with considerable variation of the ratios 
and directions of principal stresses during loading, the strains 
experimentally obtained were in good agreement with the strains pre- 
dicted by use of the deformation theory, as shown in references 9 
(p. 199), 10 (pp. 213-15) , and 11, and also in experimental work 
presented by the HA.CA at the Third Symposium on Plasticity at 
Brown University. Further experimental investigation is needed to 
determine the extent to which the variations of directions and ratios 
of principal stresses are permissible. If the variation is small, 
however, the deformation theory can be expected to give good results. 

The simplicity of the stress-strain relation given by the 
deformation theory minimizes the difficulties of solving concrete 
problems; it is therefore of both theoretical and practical importance 
to investigate the circumstances (stress-strain curves of materials, 
loading conditions, boundary conditions, and geometry of the members) 
under which the deformation theory applies. 


SIMBOLS 

The following symbols are used in this report: 

A,B,C, coefficients of nonlinear differential equations; functions 
D,E,F of a and y 

a inner radius of hole 

b outer radius of flat ring 

c outer radius of plate, very large compared with radius a 

H,J,L functions of a 

h instantaneous thickness of plate 

hinit initial thickness of plate 

l instantaneous length 
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m parameter relating to strain hardening 

n parameter relating to criterion of applicability of 

deformation theory 

r radial coordinate 

u radial displacement 

a parameter indicating ratio of principal stresses 

7 octahedral shear strain 

e logarithmic strain (or natural strain), logarithm of instan- 

taneous length divided by initial length of element 

9 angular coordinate 

O normal true stress, force per unit instantaneous area 

t octahedral shear stress 

Subscripts: 

o at inner radius of hole 

at radii b and c, respectively 
principal direction in general 

principal directions in cylindrical coordinate system 
BASIC EQUATIONS 

A thin infinite plate uniformly stressed in its plane in all 
directions and having a circular hole is shown in figure 1. The 
whole system is equivalent to a very large circular plate of 
radius c with a small concentric circular hole subjected radially 
to the same uniform stresses ■ O on the outer boundary. The 
solution obtained in such a plate within any radius b can be also 
considered as a solution of a flat ring with inner radius a and 
outer radius b loaded uniformly at the outer boundary with the 
radial stress O - 5 obtained in the plate solution. A small element 
defined by A0 and A(r+u) taken at radius (r+u) in the deformed 
state is also shown in figure 1. In the undeformed state, the, same 


b,c 

-j 

r, 0 ,z 
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element would be at radius r and defined "by &9 and Ar . Also 
shown in figure 1 are the instantaneous thickness of the element 
and the stresses acting on the element. 

The stress and strain relations based on the deformation theory 
for the case of plane plastic stress are 



where t( 7 ) is an experimentally determined function. The con- 
stants 1/2 and l/3 in equations (4) are determined from the condition 
defined by equation (l) and one of the equation (2). Only five 
of these given equations are therefore independent. 


The finite-strain concept (references 3 and 12), which considers 
the instantaneous dimensions of the element, is used because large 
deformations in the strain-hardening range will be considered. The 
stress is then equal to the force divided by the instantaneous area, 
and the strains are defined by the following equations: 
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where 2j is the instantaneous length of an element having the 
original length of ( 2 ) q and j is any principal direction. The 
octahedral shear stress-strain relation, the value of octahedral 
shear strain, and the value of the principal strains are defined "by 
the initial and final states for the paths along which the ratios of 
principal stresses remain constant during loading (reference 2)j 
thus. 




logc 


<*J>o 


or e 


e J 


l Jj 

<Vo 


The strain given "by this equation has "been referred to as "logarithmic 
strain or natural strain" (references 3 and 12) » The strain- 
displacement relations of the element shown in figure 1 are 


e 


„ d(r+u) 
dr 


(5a) 


r+U 

r 


(5h) 


e C z 


Hnit 


(5c) 


The equations of equilibrium and compatibility after the elim- 
ination of the displacements and C z (obtained in reference 2) 
are as follows: 

- € 


r 

di 


da. 




- a. 


r\ ^( € r + _ tm „ \ 

r ) r T7t \ — * (°e " °r' ® 


0 


d£ 


(6a) 



Equations (2b), (3), (4a) (4b), (6a), and (6b) are six equations 
with six unknowns, O r , C@, e r > £0> T, and 7. By using the 
transformations of references 1 (p. 189), 2, and 13, 


Qq + O r = 3 ^2* T sin a 
qq - o r ■- T cos a 
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or 

- # T 
=/\/| T 

and 

€ r = — 2L 
2V2 

£ a * 2L 

2hfZ 

The parameter a is closely related to the ratio of principal 
stresses From equation (7), 

°r _ a/T sin a - cos a 
/^3~ sin a + cos a 

This relation is plotted in figure 2 for the range of a encountered 
in the present problem. It is seen that 0 r /o g varies almost 
linear ily with a. The percentage variation of a during loading 
will then indicate directly the same percentage variation of O r /<J q. 

When this transformation is used, the equations of the problem are 
then reduced to two nonlinear differential equations with an exper- 
imentally determined relation between T and 7 . 



^sin a - a/3* cos 
^sin a + cos 0 ^ 


where 
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A 


B 


= ||^3* cos a + sin 3* sin a - cos a^ ^ 

^^3* sin a - cos ^2 

(•# 


cLT _ 7 sin a ] 1 


d 7 
7 cos 




•) 


C = 2(cos a) exp 

D = (VT sin a - cos a ) y 
E = - ^/^3~cos a + sin a) 

F = 2V2*jl - exp 7 eos a^j 

From equations (9), there is obtained 


r\ da CE - FB 

sc ■ •••• 


a/ d /r\ AE - DB 


(I) 


dl 


dy _ FA - CD 
r\ EA - BD 




J 


If a is considered as an independent variable , the calculation 
can be greatly reduced because many terms in the equations are 
trigonometric functions of a; that is, 



AE - DB 
CE - BF 


d 2, = FA - CD 
da CE — FB 


Using equations (9a) and expanding 
series gives 



( 10 ) 


into a 
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CE - BF = - 2KL - 2 tf5 EJ g(a,y ,1 g) 7)^ 



where 

H m coa a 

J = *{T sin a - cos a 


cos a + sin a 


and 


(cos a) 


?! (a-,7) = |i - e ^ 

(cos a)y 

= 1 - (cos a) 7 + I (cos 2 a )y 2 . 

S L 7 ,1 *lLz vr -Afli y 

\ * °T ZyJ T d7 V 2 J 


(lOa) 


LINEARIZATION OF EQUATIONS 


Substituting equations (10a) into equations (10) yields the 
following equations: 



-(f) 2 

7 dT 

T d7 " 

(*!z j r 0 ) r ° 



7 dT - 1 7 \y 

pa? J ^/oj 


I^Ey +± H 2 7 2 - . . .) 

2 4 / 


(lla) 
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d loge7 
da 


(-4 


1 - 


8^2 


L7 + 


LH7 2 - 


16 


1 + 


§ - (# HhK 1 - l# H? ♦ i hV - • • ) 


(lit) 

always occurs in 


It is seen in the previous equations that 

the combination 2 — . When a f (7) curve of a material is plotted 

in logarithmic scale, the slope of the curve at any value of 7 

d log e T _ 7 dT 
“ T d7 


is 


d log e 7 


The slope 


in the logarithmic plotting is a function of 7. 


7 dT 
T d7 

The use of the power function relating to the octahedral shear stress 
and strain, which approximates the actual octahedral shear stress- 
strain relation of a given material over the Btrain- hardening range, 
will be subsequently shown to give results approximately the same as 
those obtained by using the actual octahedral shear stress-strain 
curve of the material. The power-function relation (or power- law 
approximation) can be represented by the following equation: 

T (7) * Kt® 


or 


2 

T d7 


a: m 


(12) 


where E and m are constants through the strain-hardening range . 
If equation (12) is wed only to represent the T and 7 relation 
for the range of 7 along the radius of the plate, accurate results 
can be expected. It should be emphasized in such a case that the 
value of m for one material is different for different loads. 

When the value of 7 0 Is decided for any calculation, the range 
of 7 can be estimated because the values 7o/7b are Q^ite close 
for different cases (reference 2). The value of m Is then simply 
the slope of a straight line that approximates the 1(7) curve of 
the actual material plotted in the logarithmic scale for the 


range 7*, to 7 0 . 
replaced by m. 


The term 2 £1 
T dy 


in equations (11) can then be 
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A comparison of the terms in equations (11a) and (lib) with 
the numerical data obtained in reference 2 shows that 

h • • •)< 

and | - 5 . Lv + LBy 2 - , , ,) are small compared with 

V 8 /\f 2 16 / 

jl + (£) 2 % 1, jl + Vs" £ z |lj, and 1, respectively. 


and 1 , respectively. 


Equations (11) can therefore be linearized by replacing the term 

kilt 1 _ 2 L by a constant Ct . whioh is so determined that the effect 
*V 2 J 7 „ 


on £ 1 ° ge i&j due to replacing i Of i _ 2 L by C and neglecting 


2 J 7 0 


the terms 


as the effect on 


i^!^ Hy + i H 2 y 2 - . , in equation (lla), as well 

it on — ?:.°3. e .-7. due to replacing llf"!" 2 JL by C 1 and 
da i ^ J 7 1 


neglecting the terms f - — 2 — Ly + V!" LH 7 2 - . . .) and 

V 8tf2 16 J 

^ - 1^| H 7 + 1 H 2 y 2 - . . .j in equation (lib), is very small. 


Equations (ll) are then written in the following form: 


* a, (a . m 1+ w (n - Ciro) 

da W 1 + ffs 1 (m - Cj7 0 ) 


a iog e y 

da 


//s’ f - C l? c ) 


(13a) 


(13b) 


n . a - C x 7 0 


(14) 
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In equations (15) all terms under the integral sign are functions 
of a except n, which is a constant. The values of* r/a and 
yh 0 corresponding to a for various values of n can therefore 

easily he calculated by numerical integration. The constant Cj in. 
equation (14) can be determined either by comparing the variations 
of a and 7 / 7 0 with r/a for different values of n calculated 

from equations (15) with those obtained in reference 2, or by com- 
paring the terms in equations (ll). . (Because the relations of the 
variables a, y/7 0 , and r/a for different cases obtained in 
reference 2 are very similar, any available relations of these 
variables for the problem can be used to compare the terms in equa- 
tions (11) in determining the constant C^.) The value of. was 

found to be 0.7 by the first method. Thus, 

n = m - 0.7 7 Q 


(14a) 
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With this method, the values of r/a and yJy Q corresponding 

to a for a number of values of n can easily be calculated and 
given in both curve and tabular form to facilitate solution of this 
problem for a given material under a given maximum strain. 


DETERMINATION OF PRINCIPAL STRESSES AND STRAINS 


After the variation of a and y/y 0 with r/a is determined 
by the manner previously described, the principal stresses and strains 
can be obtained by equations (7) and (8) with the actual T (y) curve 
of a given material. It is seen from the equations that £ r /y, 

£ gfy } 0 r /*r, and Oq/t are functions of a only; they can be cal- 
culated for different values of a and given in tabular form to 
facilitate solution of any given case. Equations (7) and (8) can 
be written as . ' . 


°rA 

’eA 


(^/"3* sin a, - cos a) ^ 
~3 sin a + cos a) 


(7a) 



£ r/7 

€ e/y 


— i— (sin a - aTIT cos 

2^2 

v (sin a + aJ~3 cos 

2*(2 



a) 


J 


The value of principal stresses and Strains can be then obtained only 
by a simple multiplication. 


The method of determining the distribution of a, y/y 0 , a r , 

Og , £ r , and £ g along r/a will be designated the linearized 
solution. 

STRESS AND STRAIN CONCENTRATION FACTORS AT HOLE 

The stress and strain concentration factors are defined as the 
ratio of tangential stress and strain at the hole to the tangential 
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stress and strain at r = c, respectively, with c/a much greater 
than 1. For the previously described loading condition, the values 
of a at boundaries are as follows: 

at r m a, 

a r 

a 0 

at r = c, 


= 0 



Or * 09 


% 


rt 

2 


From equation (7), the stress concentration factor is 

(g 9 )o 'ofr* el ° l * °° 8 . To 

(°S>c t c (^ sin £ + cos T ° 

From equation (8), the strain concentration factor is 

<««>o 7o ( Bln ! + A/Tcos f) *7„ 

i I.. I. s z . . i h i . . . ' . ; 

(£(,)c r 0 (.m | + V? cos |) r ° 


(16 a) 


(16b) 


The stress and strain concentration factors were experimentally 
determined in reference 14 for a tension panel with a circular hole. 
An approximate equation of stress concentration factor for the same 
problem was given in reference 15. 


IDEALLY PLASTIC AMD APPBOXIMATE SOLUTION 
When n = 0, equations (15) become 


£ = exp <. 


j" F PCf da - I d (log e cos a)J 


V. \J a o 
r\ a 


JL = exp | 
7 ^ 


4 da 


a o 


a/z + tan cJ 
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/ r \2 cos Oq ^f3*( a “ a c ) 

) ~ cos a 6 

y ft cos (Xq + sin ctQ A/T(a 0 - a) 

— _ — — e 

7 q q 3 cos a + sin a 


(17a) 

(17b) 


For ideally plastic material when the infinitesimal strain concept 
is employed, m = 0 and y Q is infinitesimal; n is then equal to 

zero, this being the special case of n = 0. For this case, the 
same relation between r/a and a as in equation (17a) was obtained 
by Nadai (reference 1). 

The variation of a and y/y 0 with r/a for a material with 

n / 0 can be determined approximately by equations (17). This approx- 
imation, however, is used only to determine strains € r and e g 

because the strains are functions of a and y only. The stresses, 
which are functions of a and T , can be determined by the approx- 
imate values of a and y/y with the actual T (y) curve of the 

material. This method will be designated the approximate solution. 


CALCULATIONS , RESULTS, AND DISCUSSION 


In order to observe the degree of approximation resulting from 
the use of a power function representing the octahedral shear stress- 
strain relation over the whole strain- hardening range, a calculation 
is made with the following relation between T and 7 ; , 

T ( 7 ) = 126,000 y 0 ' 25 

The constants 126,000 and 0.25 are chosen in order to approximate 
the T ( 7 ) curve of Inconel X over the whole strain- hardening range 
(fig. 3). The variations of a, 7 , a r , Og, € r , and with r/a 

for 7 = 0.3 obtained by the power- law approximation are compared 

® •* . . 
with the values obtained from the actual T( 7 ) curve of Inconel X 

in figure 4 . These curves are seen to agree well. This result 

indicates that the method used to linearize Z — in equations (ll) 

T d 7 
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should introduce very little error. Also, if a simple analytical 
function of 'T(y) is desired in the analysis, this simple form 
T (7) = K^ m representing the t and y relation of a given material 
will give a very good approximation. 

The values of r/a and y/y 0 corresponding to a are cal- 
culated for n * -0.15, -0.1, 0, 0.1, 0.2, 0.3, and 0.4 from the 
linearized equations (15). The variations of a and y/y Q with 
r/a for different values of n are plotted in figure 5 and are 
tabulated in table I. The curves of a and y/y Q against r/a 
obtained in reference 2 for Inconel X and 16-25-6 are plotted in 
figure 6; the values of n for these cases are calculated from 
equation (14a) and are also indicated on these curves. The curves 
obtained from the linearized solution (fig. 5), which have values 
of n close to the values of n of the curves obtained from ref- 
erence 2, are also plotted on figure 6 for comparison. It can be 
seen that the results for the same n of the two solutions' agree 
very well for most cases. For simplicity, the solutions obtained 
in reference 2 based on the deformation theory are designated the 
exact solution. The simple relation n = m - 0.7 y Q can then be 
used as a good approximate criterion to find the distributions of 
a and y/y 0 with r/a for different materials and different max- 
imum strains. This generality has a very important consequence, 
because it leads to a general criterion of applicability of the 
deformation theory of plasticity for any material in the strain- 
hardening range for this problem. The criterion now is that if the 
value of n for a given material is constant or approximately 
constant in the strain- hardening range, the deformation theory can 
be applied to the problem for this material. For the special case 
of infinitesimal strain, the condition that n is constant reduces 
to the condition that m is constant, which is the same condition 
obtained by Ilyushin (reference l6) . 

The maximum variation of a at any r/a is less than 10 per- 
cent for a variation of n from -0.15 to 0.4 (fig. 5(a)). Because 
for most materials m increases with y, it can be seen from equa- 
tion (14a) that m and y Q affect n in an opposite manner, so 
that the variation of n with strain is small. The value of n 
for most materials and for any strain varies from -0.1 to 0.25 
(references 17 to 19); only a small part of this variation is due 
to variation with strain for a given material. The variations of a 
during loading will then be expected to be very small for tost 
materials; consequently, the deformation theory of plasticity is 
applicable within engineering accuracy to this problem. 
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It is also interesting to note the variations of 7 / 7 - with 
r/a for different values of n. From figure 5(h) the curves of 
7 ho against r/a near the hole, which is the region of interest, 
are close for values of n from -0.15 to 0 . 4 . 

The values of 0 r A > 0 gfa > € r f j , and Zq/j were calculated 
by equations (7a) and ( 8 a) for the range of a considered 
(table II) . The distributions of principal stresses and strains 
were obtained by the linearized solution for the following cases: 

7 0 = 0.1871, Inconel Xj y Q = 0.3000, Inconel X; and y Q = 0.4600, 
16-25-6. The results are plotted in figure 7 with results obtained 
in reference 2 for comparison. The examples y 0 = 0.3000, Inconel X, 
and 7 0 = 0.4600, 16-25-6 are chosen because they show the greatest 
deviation in a and yho between, the linearized and the exact 
solutions. The example of y Q = 0.1871, Inconel X is chosen because 
it has an average agreement between the two solutions. The 7 value 
of n for this case is equal to 0.151; a and 7/70 at cor- 
responding values of r/a can be obtained from table I by linear " 
interpolation. The distributions of principal stresses and strains 
obtained by the linearized solution agree very well with those 
obtained by the exact solution even for the case of 7 0 = 0.4600, 

16-25-6, for which the agreement in y/y 0 between the two solutions 
(fig. 6 ) is least. 

Calculations were also made for the case of ideally plastic 
material (T = constant) with the infinitesimal strain concept. 

The variations of a and 7/7 with r/a are the same as for the 

case of n = 0 (fig. 5). The distributions of principal strains 
were calculated for y Q equal to 0.1871 and 0.4600 (fig. 8). The 
curves obtained by the exact solution for y Q = 0.1871, Inconel X, 
and 7 q = 0.4600, 16-25-6 are also plotted in this figure. These 
two cases are chosen because the a and 7/70 curves deviate most 
from the curves of n = 0 (fig. 5)„ From the curves of a and t/t q 

against r/a of n = 0 and the actual T(y) curve of Inconel X and 
16-25-6, the distributions of principal stresses are calculated for 
the cases of = 0.1871 and f 0 = 0.4600, respectively; this 
procedure gives the approximate solution. 

The principal stresses are also calculated for the ideally 
plastic materials. In this case T = constant and is chosen equal 
to the value of t q for which the solutions are being compared. 

These distributions of principal stress along the radius, as well 
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as those obtained by linearized and exact solutions are plotted for 
comparison in figures 9(a) and 9(b) for y Q = 0.4600, 16-25-6, and 
for 7 0 = 0.1871, Inconel X, respectively. It can be seen from 
figure 8 that the principal strains obtained by the ideally plastic 
material with the infinitesimal strain concept (which is the same 
case as the approximate solution) gives good approximate results. 

As shown in figure 9, however, the principal 'stresses, especially 
the tangential stress (the most important part of the solution), 
obtained by the ideally plastic material are not good enough to 
approximate the actual results no matter which value of T is 
used. With slight modification as described in the approximate 
solution (using a and 7/y 0 for n = 0 and the actual T (7) 
relation of the material) , the principal stresses, especially the 
tangential stress, obtained agree much better with the exact 
solution. This simple approximate method can therefore give good 
results. It also indicates that the ratio of the strains along 
the radius to the maximum value and the ratio of principal stresses 
are essentially independent of the T (7) curve of the material and 
the maximum strain of the plate, but the stresses are very much 
dependent on the t( 7) curve of the material. This result not only 
further confirms the conclusions obtained in reference 2, but also 
generalizes those results to most materials (that is, not limited by 
the 1(7) curves of Inconel X and 16-25-6 previously investigated). 

All the results and discussions are true only for this problem 
under plastic deformation in the strain-hardening range in which the 
elastic strains are small compared with the plastic strains, and 
neither time and temperature effects nor unloading are present. 


CONCLUSIONS 

The results obtained for a thin plate with circular hole in the 
strain- hardening range, in which the elastic strain is small compared 
with the plastic strain, show that: 

1. The results obtained by the linearized solution agree very 
well with those obtained by the exact solution based on the deformation 
theory of plasticity. The amount of computation required is very 
much reduced by the linearization. With, the tables or curves pre- 
sented herein, the distributions of octahedral shear stress and strain 
as well as the principal stresses and strains can be obtained for 
a given material under a, given maximum strain of the plate by a simple 
multiplication. s 
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2. The variation of a parameter determined from the octahedral shear 
stress-strain curve of the material can he used as a simple general 
criterion of applicability of deformation theory for this problem. 

3. The results previously obtained for Inconel X and 16-25-6 

were extended to most materials: namely, the ratios of the strains 

along the radius to the maximum value and the ratio of principal 
stresses are essentially independent of the octahedral shear stress- 
strain curve of the material, but the stress distributions are very 
dependent on the material. 

4. The results obtained for the ideally plastic material with 
the infinitesimal strain concept give good approximate values of 
principal strains, but not of principal stresses. 

5. Sufficiently accurate values of principal stresses can be 
obtained by the approximate method of using the strains obtained by 
the ideally plastic material, together with the actual octahedral 
shear stress-strain relation of the material. 

6. If a simple analytical function representing the octahedral 
shear stress-strain relation is required for analysis, the power- law 
approximation can be used, 

Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 

Cleveland, Ohio, September 20, 1950. 
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TABLE II - VALUES OF C T fj, Gg/ 7 , O r / T, AMD O g/t 


FOR VARIOUS VALUES OF a 



a 


£ g/7 

V* 

V 

0.5236 

-0.3536 

0.7071 

0 

2.1214 

.5934 

- .3100 

.7054 

.1710 

2.2016 

.6632 

- .2649 

.7002 

.3411 

2.2713 

.7505 

- .2067 

.6890 

.5511 

2.3424 

.8378 

- .1470 

.6725 

.7571 

2.3961 

.9250 

- .0862 

.6509 

.9571 

2.4312 

1.0123 

- .0247 

.6243 

1.1502 

2.4481 

1.0821 

.0246 

.5997 

1.2980 

2.4480 

1.1519 

.0739 

.5721 

1.4397 

2.4361 

1.2130 

.1167 

.5456 

1.5582 

2.4160 

1.2741 

.1590 

.5172 

1.6705 

2.3868 

1.3265 

.1949 

.4912 

1.7622 

2.3547 

1.3788 

.2302 

.4639 

1.8486 

2.3160 

1.4137 

.2534 

.4451 

1.9036 

2.2870 

1.4486 

.2763 

.4255 

1.9562 

2.2547 

1.4661 

.2876 

.4156 

1.9818 

2.2377 

1.4835 

.2988 

.4056 

2.0065 

2.2201 

1.4923 

.3044 

.4005 

2.0186 

2.2109 

1.5010 

.3100 

.3954 

2.0309 

2.2017 

1.5097 

.3155 

.3903 

2.0426 

2.1920 

1.5184 

.3210 

.3851 

2.0543 

2.1825 

1.5272 

.3265 

.3799 

2.0659 

2.1727 

1.5359 

.3320 

.3747 

2.0774 

2.1629 

1.5446 

.3374 

.3695 

2.0887 

2.1528 

1.5533 

.3428 

.3642 

2.0996 

2.1424 

1.5621 

.3482 

.3589 

2.1107 

2.1320 

1.5708 

.3536 

.3536 

2.1214 

2.1214 
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Figure 2. - Variation of ratio of principal stresses o v /oq with par- 
ameter a. 
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(d) Variation of principal strains with proportionate 

radial distance. 


Figure 4* - Concluded. Comparison of results obtained for t(y) curve of 

0.25 

Inconel X and power-function approximation t(y) * 126,000 y * * 

Y. = 0 • 30 . 
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Proportionate radial distance, r/a 

(a) Variation of parameter o with proportionate radial distance, r/a. 

Figure 5. - Variation of parameter a and the ratio of octahedral shear strain to octahedral shear strain at hole 
v/y with proportionate radial distance r/a obtained by linearized solution for various values of n. 
o 
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Figure 5. - Concluded. Variation of parameter a and the ratio of octahedral shear strain to octahedral shear 
strain at hole y/y with proportionate radial distance r/a obtained by linearised solution for various val- 
ues of n. ° 
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(a) Variation of parameter a with proportionate radial distance; Inconel X* 

Figure 6. - Comparison of variation of parameter a and ratio of octahedral shear strain to octahedral shear 
strain at hole y/y with proportionate radial distance r/a obtained by linearized solution and by exact 
solution. 0 
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(c) Variation of ratio of octahedral shear strain to maximum octahedral shear strain at 
hole y/y 0 with proportionate radial distance; Inconel X. 

Figure 6. - Continued. Comparison of variation of parameter o and ratio of octahedral shear strain to octa- 
hedral shear strain at hole y/y with proportionate radial distance r/a obtained by linearized solution and 
by exact solution. ° 
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(d) Variation of ratio of octahedral shear strain to the octahedral shear strain at 
hole y/y 0 with proportionate radial distance; 16-25-6* 

Figure 6. - Concluded. Comparison of variation of parameter a and ratio of octahedral shear strain to octa- 
hedral shear strain at hole y/y 0 with proportionate radial distance r/a obtained by linearized solution and 
by exact solution. 
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(a) Variation of principal strains with proportionate radial distance. 

Figure 7. - Comparison of principal stresses and strains obtained by linearized solution and exact solution. 
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Proportionate radial distance, r/a 

(b) Variation of principal stresses with proportionate radial distance, 
figure 7. - Concluded. Comparison of principal stresses and strains obtained by linearized solution and exact solution. 
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Figure 8. - Comparison of distributions of principal strains along proportionate radial distance obtained by approxi- 
mate solution and ideally plastic with exact solution. 
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(b) — 0.1871 , Inconel X. » 

Figure 9. - Concluded. Comparison of distributions of principal stresses along proportionate radial 
distance obtained by different methods. 
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